function result = sigmasquarebar(m,v)

% result = quadgk(@(s) s.^2.*normpdf(s,m,v),-Inf, Inf);
% disp(m^2+v^2)
% result = quadgk(@(s) s.^4.*normpdf(s,m,v),-Inf, Inf);
% disp(m^4+6*m^2*v^2+3*v^4)
result = quadgk(@(s) (2+sin(s)).^2.*normpdf(s,m,v),-Inf, Inf);

space = linspace(-10,10,10^6);

f=@(s) (2+sin(s)).^2.*normpdf(s,m,v);

disp((sum(f(space))-f(-10)/2-f(10/2))*20/10^6)